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Abstract: We illustrate how QCD color structure elegantly can be decomposed into orthogonal multiplet 
bases corresponding to irreducible representations of SU(W) with the aid of Wigner 3j and 6j coefficients. 
We also show how to calculate the relevant 3j and 6j coefficients using multiplet bases and birdtrack 
techniques and argue that only a relatively small number of Wigner 3j and 6j coefficients are required. For 
up to six gluons plus quark-antiquark pairs we explicitly calculate all 6j coefficients required for up to NLO 
calculations. 
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1 Introduction 

One of the challenges associated with the high multiplicity of color charged particles at the LHC is the 
treatment of the non-Abelian color structure in QCD. The traditional method for mastering this issue is 
to use non-orthogonal bases, most notably trace bases [1~12], where the basis vectors are products of open 
and closed quark-lines, and color flow bases [13], where the adjoint representation of the gluon has been 
rewritten in terms of quark indices. 

While these strategies come with several advantages — the conceptual simplicity, the simplicity of 
gluon emission, the simplicity of gluon exchange [8], the natural interpretation in terms of flow of color [13] 
and the existence of recursion relations for scattering amplitudes [5, 14-19] — the non-orthogonality (and 
the overcompleteness) is a very severe drawback when it comes to squaring amplitudes for processes with 
many colored external legs. 
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Figure 1. Example of a six-gluon multiplet basis vector (a), and basis vectors for one gg-pair and five gluons (b,c). 
In the six-gluon case the two first incoming gluons are forced to be in representations ai etc. 

In order to cure this, a general recipe for the construction of minimal and orthogonal basis vectors has 
recently been proposed [20]. For these bases the color space squaring of basis vectors scales only as the 
dimension of the vector space, i.e., roughly as the factorial of the number of gluons and gg-pairs in the limit 
where the number of colors, N^, goes to infinity, and roughly as an exponential for finite Nc, rather than 
the square of the number of spanning vectors [20]. These bases thus have the potential to speed up exact 
calculations in color space very significantly. However, to realize this potential, it remains to argue that 
the decomposition of scattering amplitudes into multiplet bases can be accomplished efficiently. 

In this paper we take important steps in this direction by showing how Feynman diagrams can be 
decomposed into multiplet bases in an efficient manner using Wigner 3j and 6j coefficients, which can be 
calculated and stored once for all, see for example [21, 22]. Furthermore we argue that even for processes 
with very many external legs we need only a manageable number of 3j and 6j coefficients and that, using 
the basis from [20], these can be calculated. 

This article is organized as follows: In section 2 we recapitulate the properties of multiplet bases (as 
constructed in [20]). Following this we illustrate how the decomposition into these bases can be efficiently 
achieved using 3j and 6j coefficients in section 3. Then we show how to calculate all relevant 3j and 6j 
coefficients in section 4, and argue that the number of such coefficients is very manageable. For leading 
order processes with up to four gluons plus gg-pairs we explicitly show all needed 6j coefficients, whereas 
all 6j coefficients required for up to six gluons plus gg-pairs for leading and next to leading order processes 
are electronically attached. Finally we conclude and make an outlook in section 5. 

2 Multiplet bases 

Orthogonal bases corresponding to irreducible representations of subsets of particles may in principle be 
constructed in many different ways. In the present paper we follow the construction in [20], and start by 
considering multiplet bases for pure gluon processes. 

Multiplet bases are based on sub-grouping partons into states transforming under irreducible repre¬ 
sentations of SU(A^c)- More specifically, we may divide Ng gluons into [^] “incoming” gluons and 
“outgoing” gluons, and then subgroup the gluons on the incoming and outgoing side, and force the sets to 
transform under irreducible representations, as illustrated in figure 1(a), where double lines denote arbitrary 
representations. 
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Table 1. The dimension of the full vector space (all orders) for Nf. = 3 and in the iV^ —t oo limit. In the case of 
gluons only the dimension of the relevant space can be further reduced by imposing charge conjugation invariance. 


For the treatment of quarks, we note that for each incoming quark (outgoing antiquark) there is an 
outgoing quark (incoming antiquark). The quarks and antiquarks may thus be grouped into pairs which 
together transform either as singlets or as octets, since 3(8)3 = 1©8^. Each qq-paiv transforming under the 
singlet representation corresponds to a (5'^^function in color space, whereas a gg-pair forming an octet can 
be replaced by a gluon index using the SU(A^c) generators (t^)'^g. Knowing orthogonal bases for processes 
with an arbitrary number of gluons we can therefore construct bases for Nqq gg-pairs by letting Nqq of the 
gluons split into a gg-pair or be removed and replaced by 6‘^q. 

Being constructed to sort groups of partons into states transforming under irreducible representations, 
the multiplet basis vectors have no direct relation to the order in perturbation theory; the bases are valid 
to all orders. Instead the bases can trivially be made minimal for a given Nc by crossing out basis vectors 
corresponding to representations which only appear for higher Nc- The dimension of the full vector space is 
thereby reduced from growing roughly as a factorial in the number of gluons plus gg-pairs, to growing as an 
exponential [20], cf. table 1. While this reduction in dimension may become significant for a large number 
of colored partons, the main gain is clearly the orthogonality, reducing the number of terms appearing in 
the squaring of an amplitude from the square of the number of spanning vectors to the number of spanning 
vectors. 


^We use the terminology octets, 27-plets etc. even for Nc 7^ 3, although the dimensions of the representations clearly vary 
with Nc- The definitions of all representations relevant for this paper, for general Nc, can be found for example in section 1.3 
and 4 of [20], where the Young diagrams of the representations and the smallest Nc for which the representations exist are given. 
The general Nc dimensions of all encountered representations can be found in the attached .m-file, WignerCoeflicients.m. 


- 3 - 










Figure 2. Schematic example of scalar product between a vector in the six gluon multiplet basis and a general 
color structure, where the dashed lines indicate color contractions. 

3 Decomposition 

Due to confinement we are only interested in color summed/averaged quantities, and it is not hard to argue 
that the action of summing over all external indices corresponds to evaluating a scalar product. Letting Ci 
and C 2 denote the color structures of two amplitudes the scalar product is given by 

(C1|C2)= er^a,...^a,a,... (3^^) 

ai, a2,... 

with ai = 1, ...,Nc if parton i is a quark or antiquark and o, = 1,— 1 if parton i is a gluon. 

As scalar products correspond to color structures where all quark, antiquark and gluon indices have 
been contracted, we may graphically represent them as fully contracted vacuum bubbles, as depicted in 
figure 2. 

One way of performing the color sum goes via removal of four-gluon vertices in favor of triple-gluon 
vertices, rewriting of structure constants in terms of traces, usage of the Fierz identity (completeness 
relation) to remove gluon propagators, and finally counting of closed quark-loops, see for example [11, 12]. 
While this scalar product evaluation procedure scales much better than summing over explicit indices, the 
(potential) doubling in the number of terms with each structure constant and each Fierz identity, makes 
the evaluation of scalar products between Feynman diagrams and basis vectors relatively expensive for 
amplitudes with many (external and internal) gluons. We will see below that using the birdtrack method 
we both reduce the number of terms and avoid calculating scalar products with a large class of basis vectors 
for which the projection vanishes. 

We thus suggest a strategy based on repeated usage of the completeness relation 
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and Schur’s lemma 



for performing color contractions for arbitrary representations, see [21] for a good introduction. Using the 
completeness relation, eq. (3.2), and Schur’s lemma, eq. (3.3), on a vertex correction gives 


c 



The sum in eq. (3.4) runs over all instances of the representation 7 in a ( 8 ) (3, corresponding to different 
vertices a. For example, if a, /3 and 7 are all octets, the right hand side may contain vertices with both 
the symmetric and antisymmetric structure constant. In general for QCD (any Nc), if one of a, f3 and 7 is 
the octet representation, the sum runs over up to two {Nc — 1 ) vertices if the other two representations are 
equal. If they are different the sum only contains one term [20]. 

The goal of the strategy is to express vacuum bubbles, for example the one in figure 2, in terms of 
vacuum bubbles of the form in the numerator and denominator in eq. (3.4), known as Wigner Qj and 3j 
coefficients, respectively. A vacuum bubble will contain loops of representations; as it is fully contracted 
every path through the bubble must close. For a loop of representations with n vertices, completeness 
relations can be applied to relate the loop to a sum over loops with fewer vertices. For example, considering 
a loop with six vertices, one possibility of applying the completeness relation is 



where we suppress representation and vertex indices. This step of the procedure is independent of the 
number of vertices in the loop on the left hand side of eq. (3.5). Hence it can be applied repeatedly to 
an n-vertex loop until the expression is a sum over loops with three vertices. The remaining three-vertex 
loop, i.e., vertex correction, can then be removed using eq. (3.4). This step removes two vertices, since a 
three-vertex loop is replaced by a single vertex. As this procedure is independent of the number of vertices 
in the loop, loops can be repeatedly contracted until any vacuum bubble equals sums over Wigner 3j and 
6 j coefficients. 
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3.1 An explicit example 

As an example, let us denote gluons by plain lines without arrows and consider the color structure of the 
Feynman diagram 



(3.6) 


for qq —?■ bg. From a color structure perspective, any leg can be changed from the incoming side to the 
outgoing and vice versa. Hence eq. (3.6) may alternatively be drawn as 




(3.7) 


In eq. (3.7) the gg-pair is in the adjoint representation, therefore it is orthogonal to all basis vectors with 
the (^g-pair in the singlet representation (cf. figure 1) by the tracelessness of the generators. 

The vacuum bubble that has to be evaluated to express eq. (3.6) in a multiplet basis is 



where the equality comes from applying Schur’s lemma to remove the two-vertex loop with the quarks. The 
smallest loop on the right hand side of eq. (3.8) is the vertex correction involving the representation ai, it 
can be removed using eq. (3.4), resulting in 



The shortest loop in the large vacuum bubble on the right hand side has four vertices and there are several 
possible choices of four-vertex loops, one is 



(3.10) 
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Contracting this loop gives 


b c 



Replacing the four-vertex loop in eq. (3.9) with the right hand side of eq. (3.11) results in 


(3.11) 



(3.12) 


The final expression is, by drawing the rightmost vacuum bubble in the same shape as the other 6j coeffi¬ 
cients, 



While the sum over vertices in eq. (3.14) in principle may contain more than one term, in all cases en¬ 
countered here, the vertices have been chosen such that only one term contributes, giving at most a minus 
sign. 

From eq. (3.13) one can immediately determine that this color structure is orthogonal to many of the 
basis vectors. It is orthogonal to all of the basis vectors where the qq-paii is in a singlet, as mentioned. It 
is also orthogonal to all basis vectors where 02 0 8 (8) 8, as seen from the vertex labeled a on the right-hand 
side of eq. (3.13). Since the orthogonality is manifest from constraints on the representations, the projection 
onto many basis vectors need not be calculated, which we expect to significantly speed up computations. 
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The right hand side of eq. (3.13) can now be explicitly evaluated. As an example the vacuum bubble 
for the (unnormalized) basis vector with the representations ai = 27, 02 = 8 and 03 = 1 is 



where the values of the Wigner 3j coefficients have been normalized to one, c.f., eq. (4.13) and eq. (4.14), 
and the 6j coefficients are taken from table 3. Using the standard normalization of the antisymmetric 
triple-gluon vertices in the original color structure, we should multiply with a factor (y^ 2 A^c 2 r(A"| ~ 1 ))^) 
where Tr is defined by the generator normalization, Tr[t“t^] = and typically chosen to be 1/2 

or 1. Similarly, the 3j coefficient from the vertices would in the standard normalization give a factor 
Tr{N^ — 1), rather than 1. If we want to project onto normalized basis vectors we also need to multiply with 
the square roots of the dimensions of the representations ai, q ;2 and 0:3 and a factor I/^/Tr from the qq- 
pair giving an additional factor \f\\/{N'^ — 1)^+ 2Nc — 3)/4/-v/T/j. We note that the quadruple 
sum in eq. (3.13) typically only contains a few non-zero terms, making the above procedure for scalar 
product evaluation by far superior to the method of rewriting all representations in the fundamental and 
anti-fundamental representation, currently employed in for example [9, 11, 12]. 

We also remark that this procedure of decomposing color structure is entirely general. In the following 
we focus on leading order (LO) and next to leading order (NLO) processes, but there is no conceptual 
difference in treating color structure at higher orders. 

An increasingly popular strategy for evaluating amplitudes is to use recursion relations of various kinds 
rather than plain Feynman diagrams. While the present paper addresses generic (Feynman diagram) color 
structure decomposition, it is also of interest to study how efficient multiplet basis color decomposition can 
be achieved in recursive approaches. This is done in detail in [24], for the case of BCFW [18, 25] recursion 
of tree-level MHV gluon amplitudes, where it is argued that the recursion can be formulated and efficiently 
implemented directly in the multiplet basis (cf. table (1) in [24]). Using this recursion, the amplitudes can 
— if desired — also be expressed in terms of color-ordered amplitudes (see for example eqs. (2.15) and 
(B.5) in [24]). 


4 Wigner 3j and 6j coefficients 


The procedure of section 3 expresses any vacuum bubble as a sum over factors of Wigner 3j and 6 j 
coefficients. The Wigner 3j coefficients are easily evaluated by contracting two vertices. 



(4.1) 


By construction the Wigner 3j coefficients of two different instances of the vertex between three represen¬ 
tations must be zero in order for the completeness relation to hold. 
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In the initial vacuum bubble there are only vertices where at least one representation is the adjoint 
representation. As is argued in appendix A, the loops for the contraction strategy of section 3 can be chosen 
such that the completeness relations are only applied where at least one of /x and v in the completeness 
relation eq. (3.2) is the adjoint representation. This results in that every new vertex introduced from the 
completeness relation also contains at least one adjoint representation. There are only two different types of 
Wigner Qj coefficients where every vertex has at least one adjoint representation that cannot be simplified 
to Wigner 3j coefficients by applying Schur’s lemma, 



0 0 


and 

(4.3) 

0 

In eq. (4.2) the black blob in the left Qj coefficient is the antisymmetric structure constant of SU{Nc), 
j^jabc _ ^ I'rjij, _ 'j'j. ^ j^j^g white blob is the symmetric structure constant, = 

^3^ [Tr + Tr . 

Internal quark-lines can be handled by applying the completeness relation and rewriting traces over 
three generators in terms of and A quark loop with an arbitrary number of gluons attached can 
be rewritten as 




has been used. By applying this repeatedly, an arbitrarily large quark loop can be expressed only in terms 
of structure constants. This introduces vertices, the last term in eq. (4.4), which in turn require the 
second coefficient in eq. (4.2). 
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4.1 Evaluation of Wigner 6j coefficients 


Given a set of basis vectors from the basis construction described in [20], the Wigner Qj coefficients can be 
calculated. The Wigner 6j coefficients of the form of eq. (4.3) can be evaluated by contracting two basis 
vectors which differ only by one representation and where two gluons are crossed, 



Here Schur’s lemma can be applied repeatedly on the two-vertex loop in the middle, and on the two-vertex 
loop formed by the trace, giving factors of Wigner 3j coefficients over dimensions of representations. Hence 
eq. (4.6) is equal to 


“r5i-2 




(4.7) 


where uq an -2 are taken to be the octet representation. To note is that if the representations 
and appear in basis vectors with fewer gluons, then these basis vectors can be used instead. 

For the calculation of Wigner 6j coefficients of the form of eq. (4.2), for n = Ng + Ngq partons, n-gluon 
basis vectors are contracted with (n — l)-gluon basis vectors. The difference of one external gluon between 
the two vectors is required since three gluons are to be contracted with one or 4“^'^ vertex. In the 
following a grey blob is used as a placeholder for the and the vertices. The vectors are contracted 
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as 



where oq and a „_3 are the octet representation. That it is sufficient with an (n — l)-gluon and an n- 
gluon basis vectors is due to constraints that can be placed on the required Wigner coefficients, derived in 
appendix A. If the Feynman diagram to be decomposed contains quarks, the grey blob is either or 

^abc, ■£ 'I- contains gluons the Wigner coefficients with are not required. 

One might worry that the value of the Wigner coefficients of the higher representations could depend 
on the construction history of the vertices, such that, for example, the Wigner 6j coefficients containing 
a vertex with a 35-plet with the construction history a = (10,35) would differ from those with vertices 
constructed from a = (27,35). In the case of unique vertices, i.e., vertices between three representations 
which can only be combined to a vertex in one way, Schur’s lemma and isomorphism guarantees that this 
can not happen (as long as the vertices have the same normalization and matching sign conventions). In the 
case of vertices that appear in several instances, for example the two vertices involving the representations 
35, 35 and 8, one can prove that these sets of vertices can be chosen to give identical 6j coefficients. In 
appendix B we give the birdtrack proof of both of these statements. 


4.2 Required Wigner coefficients for LO and NLO color structures 

The Wigner coefficients of eq. (4.2) and eq. (4.3) have many symmetries which can be used to reduce the 
number of coefficients that has to be calculated. The symmetries of the relevant Wigner 6j coefficients, 
proven in appendix C, are: 
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(i) Rotation symmetry 



where the * by the vertices is there to indicate that each vertex is to be understood as the conjugated 
version of the vertex on the left hand side. Using the completeness relation we note that the conjugated 
vertex may in principle be a linear combination of various vertices, 


^ a P 

V 



(4.11) 


Note that even for real representations, conjugating internal lines in the vertex may change it. All 
vertices encountered in the present paper have been chosen such that the sum above only contains one 
term, meaning that conjugating gives at most a minus sign, as for the triple-gluon vertex. In general 
this also holds whenever a /3, since in a ( 8 > 8 , there is at most one instance of the representation /3 
[ 20 ]. 

(iii) Reversion symmetry 



With the symmetries stated in eqs. (4.9-4.12), along with the constraints on the representations a, 
P, 7 and 5 proven in appendix A, an upper limit for the required number of Wigner coefficients has been 
calculated, and is shown in table 2. We note that, in comparison to the number of basis vectors (see table 1), 
the number of needed 6 j coefficients — which can be calculated once and for all — grows very slowly with 
the number of partons. It should nevertheless be remarked that the numbers in table 2 necessarily are 
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n — Ng + Nqq 

4 

5 

6 

7 

8 

9 

10 

11 

12 

LO gluons Nc = 3 

21 

39 

106 

152 

254 

318 

452 

536 

705 

NLO Nc = 3 

29 

55 

120 

176 

272 

350 

476 

576 

733 

LO gluons Nc ^ 00 

28 

68 

313 

636 

1 777 

3 095 

7 289 

12 009 

25 487 

NLO Nc^oo 

44 

108 

389 

808 

2 023 

3 693 

8 077 

13 783 

27 613 


Table 2. Upper limits for the needed number of 6j coefficients for tree-level gluon amplitudes and NLO color 
structures with up to n gluons plus gg-pairs. These numbers are calculated using the symmetries in eqs. (4.9-4.12) 
and the representation constraints from appendix A. The actual numbers of non-vanishing 6j coefficients depend on 
the choice of (generalized) vertices. 


upper estimates since the 6j coefficients depend on the choice of generalized vertices. For the case when a 
representation 7 appears more than once in the tensor product of two representations a®/?, such that there 
are several vertices connecting the representations a, (3 and 7 , there is in principle an infinite number of 
ways of defining the corresponding vertices. Since the 6 j coefficients depend on the vertices, this is reflected 
in the values of the 6 j coefficients. In particular, a clever choice of vertices may lead to the vanishing of 
some 6 j coefficients which would not vanish for another choice. 

Before actually evaluating any 6 j coefficient we need to decide on the generalized vertices connecting 
three general representations. We choose our convention such that all vertices are normalized to one in the 
sense that the corresponding 3j coefficient is one. 



= 1 (for values of 6 j coefficients). 


(4.13) 


Especially we remark that this normalization is applied to the antisymmetric triple-gluon vertex 


giving the vertex 


e- (for values of 6j coefficients). 


(4.14) 



1 

^2N,Tr{N^-^ 



(for values of 6j coefficients). 


(4.15) 
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Table 3. Non-vanishing Wigner 6j coefficients required for up to NLO QCD color structures with up to four 
external gluons plus qq-pairs. Note that all vertices (in particular the antisymmetric triple-gluon vertices) have been 
normalized such that the corresponding 3j coefficient is one, c.f. eq. (4.13) and eq. (4.14). 
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As an example, the non-vanishing Wigner 6j coefficients required for QCD color structures appearing 
in calculations up to NLO with up to four external gluons plus g'g-pairs are shown in table 3. Vanishing 6j 
coefficients, 6j coefficients related by symmetries, and 6j coefficients only required for Nc > 4: have been 
omitted. The sign conventions for the vertices involving two different representations of 10, 10, 27 and 0 
and one octet are given by normalizing the vertices 


a 



oc 


fi 


(4.16) 


s.t. the normalization constant is positive (and fixed by requiring the 3j coefficient to be one). 

For the two vertices involving 27, 8, and 27 one vertex is chosen antisymmetric (a) under exchange 
of vertex order, and is constructed as in eq. (4.16), whereas the other vertex is chosen symmetric (s), and 
proportional to 











(4.17) 


with positive proportionality constant. For higher Vc, the vertices between 0, 8, and 0, are defined similarly. 
For the decuplet vertices involving 10, 8 and 10, one vertex is chosen as in eq. (4.16), whereas the remaining 
vertex — which only appears for Vc > 4 — is given by the orthogonal linear combination of eq. (4.16) and 
eq. (4.17) with positive constant in front of eq. (4.17) and negative in front of eq. (4.16). The corresponding 
vertices with anti-decuplets, involving 10, 8 and 10, are defined as the conjugate of the decuplet versions, 
which implies a minus sign in front of eq. (4.16), for both cases. These signs are also encoded in the 
electronically appended five-gluon basis. 

The sign conventions of the other vertices, which appear twice in the 6j coefficients in table 3, do not 
change the sign of the coefficients. 

Going beyond four gluons plus gg-pairs, the full set of 6j coefficients required for up to NLO calculations 
with up to six gluons plus gg-pairs are attached in a human and Mathematica readable .m-file. Here the 
relevant sign conventions are defined by the appended six- and seven-gluon basis vectors. Since we view 
the basis vectors as constructed from vertices, the usage of the birdtrack method imposes sign correlations 
between the basis vectors. For this reason — and since we require that conjugating a vertex should give the 
vertex with conjugated representations when possible — some six-gluon basis vectors have changed sign 
w.r.t. the basis published along with [20]. 


5 Conclusion and outlook 

We have demonstrated how QCD color structure elegantly can be decomposed into multiplet bases with 
the aid of Wigner 3j and 6j coefficients, which we have shown how to calculate using the multiplet bases 
from [20]. 
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We have also argued that only a relatively small set of such coefficients are needed, and that the number 
of required coefficients is severely reduced for Nc = 3 compared to the limit W oo. For leading and 
next to leading order processes and up to six gluons plus gg-pairs we have explicitly evaluated all necessary 
coefficients. 

We remark that although the discussion in the present paper has focused on the decomposition of color 
structure associated with Feynman diagrams, the same principle can be applied in other contexts, such 
as the color structure of gluon emission for parton showers, the calculation of soft anomalous dimension 
matrices, or recursive approaches to scattering amplitudes [24], 
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A Representations in 6j coefficients 

This appendix contains a proof that Wigner coefficients of the form of eq. (4.2) and eq. (4.3) are sufficient 
for decomposing QCD color structures by the method of section 3. For the proof, specific choices of loops 
in the contraction procedure are examined. It is shown that such loops can always be found in the original 
vacuum bubble, and in every subsequent step. The last part of this appendix puts constraints on the 
representations occurring in the Wigner coefficients, for both tree-level and NLO color structures. We 
stress that the contraction strategy of this appendix should be viewed as a proof of the constraints that 
can be put on the Wigner coefficients, and not as a suggestion for implementation of the method. In actual 
examples, one can often, as in section 3.1, choose contractions such that most loops only involve two or 
three vertices. 

For this appendix it will be important to differentiate between the adjoint representations appearing 
in the color structure to be decomposed, and the arbitrary representations in the basis vector or from 
completeness relations (which may well be adjoint representations). In this appendix any mention of the 
adjoint representation will exclusively refer to the adjoint representations from the initial color structure, 
i.e., they are adjoint representations for all of the basis vectors. The labeled representations, e.g., ai, 02 
and as in figure 1, are referred to as the arbitrary representations. 

A.l Gluon-only color structures 

To completely contract vacuum bubbles with both external quarks and gluons, three different types of loops 
will be required. In the following, tree-level and NLO gluon-only color structures are handled first. The 
additional type of loops required for an arbitrary number of quarks is dealt with afterwards. The vacuum 
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bubble for LO gluon-only color structures will always contain at least two loops of the form 



where the characterizing feature is that there is only one vertex from the initial color structure (the gray 
blob, in the gluon-only case this is always The total number of vertices in the loop does not matter 

for the principle. To reduce such a loop, the completeness relation, eq. (3.2), and the contraction of a vertex 
correction, eq. (3.4), can be applied to the two red representations in 



(A.2) 


Repeating this procedure will result in a vertex correction containing the gray blob (for the loop in the 
above example, this step would need to be repeated two more times). The vertex correction with the gray 
blob gives 



E 

a 



(A.3) 


for some representations a, j3 and 7 . This last step removes two vertices, one gray blob, i.e., a vertex 
from the initial color structure and one vertex between arbitrary representations. Since every loop that is 
contracted removes one vertex of each kind, the resulting vacuum bubble is topologically equivalent to a 
vacuum bubble for a tree-level color structure with n —1 external gluons. After a loop of the form of eq. (A.l) 
has been contracted, there must thus exist at least two loops of the type in eq. (A.l) in the resulting color 
structure by the above argument. Hence any LO gluon-only color structure can be completely contracted 
by repeatedly contracting loops of the form of eq. (A.l). 

Only choosing loops of the form of eq. (A.l) is sufficient for tree-level gluon-only color structures. For 
higher orders, it is not always possible to choose loops of this form. At NLO this happens for diagrams 
where all external gluons are attached to the loop, such as 
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(A.4) 















For all other gluon-only NLO color structures there is at least one with two uncontracted indices, 
meaning that a loop of the from in eq. (A.l) can be found, such that it is possible to contract loops as in 
eq. (A.2). For color structures of the form of eq. (A.4), there always exists loops of the form 



(A.5) 


Similarly to the loop in eq. (A.l), the steps detailed in eq. (A.2) remain valid. However, at the end, instead 
of contracting a loop of the form in eq. (A.3), a loop with four vertices is encountered. 



(A.6) 


Handling a loop like this removes two vertices from the color structure, and none from the chain of arbitrary 
representations. Since an NLO color structure has two vertices more than a tree-level color structure, the 
topology after the contraction (remembering that the two triple-gluon vertices of eq. (A.5) are in the loop 
of the NLO color structure) is equivalent to that of a tree-level color structure. Note that a loop of the type 
in eq. (A.5), need not be the first loop to be contracted (most NLO diagrams contain no loop of the form in 
eq. (A. 4)), but might at some step be encountered, and necessary to contract to continue the contraction 
of the vacuum bubble. In fact, all loop contractions apart from one, where eq. (A.5) is inserted to break 
the NLO loop, can be of loops of the form in eq. (A.l). In this way, we can thus contract any gluon-only 
NLO diagram. 

Color structures of arbitrary order in perturbation theory can be decomposed by contracting similar 
loops to eq. (A.l) and eq. (A.5), in general with three or more vertices from the color structure, instead 
of one or two, respectively. These have been avoided above in order to be able to put constraints on the 
representations appearing in the Wigner coefficients. 

When dealing with internal quark loops, the terms coming from the first term on the right hand side 
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of eq. (4.4), appearing when removing the internal quark loop, will give rise to loops of the form 



(A.7) 


Also for this loop, the first step for the first loop type, eq. (A. 2), can be performed analogously. The 
difference again occurs in the last step, where a vertex correction of a different type remains of the loop. 



(A.8) 


In this step two vertices from the arbitrary representations are removed, and none from the initial color 
structure which, however, already, for these “singlet” terms, contain two vertices less. Thus contracting loops 
of this type allows for gluon-only color structures with an internal quark loop to be completely contracted 
as well. 


A.2 Color structures with external quarks 


For color structures with quark-antiquark pairs, the color structure after using eq. (4.4) is not guaranteed 
to have any loops of the form of eq. (A.l) or eq. (A.7). As will be argued below, the consequence of this 
is that loops of the type in eq. (A.5) are required also at LO. When the quarks are combined into octet 
or singlet representations in the basis vectors to match the basis vectors in figure 1 (b) and (c), they will 
form traces over some number of generators. These quark traces can be simplified by using completeness 
relations and removing loops. As we will see, after the quarks have been removed, an amplitude containing 
only gluons and and vertices remains. This gluon-only color structure will be of the same order 
as the starting color structure with quarks (or lower), hence it can be contracted as described in section A.l 
for LO and NLO color structures. 


The color structure of a general QCD amplitude with Ng gluons and Ngq quark-antiquark pairs will be 
of the form 



19 









where the gray blob can be of arbitrary order in perturbation theory. Considering one of the external 
quarks, with some label the amplitude can be written as 

f*' * 


j 

where j is one of the external antiquarks and the dots to the right represent that there can be any number 
of gluons connecting to the rest of the color structure. 

To manipulate eq. (A. 10), a completeness relation, eq. (3.2), can be applied to the open quark-line with 
i and j (analogously to eq. (4.4)) giving 



(A.IO) 





(A.ll) 


where the big blob has been suppressed and the small gray blob represents + Repeated application 
of eq. (A.ll) allows for the color structure to be written as a sum over different contributions where quark 
i and antiquark j are in either a singlet or an octet, i.e. eq. (A.IO) can be written as 



(A.12) 


where the thin double line a is either a singlet or an octet. If this step is performed for every incoming 
quark, the total color structure is of the form 



The order of the gray blob is the same as for the original gray blob in eq. (A.9), or lower, since eq. (A.ll) 
does not introduce any loops, but contains a singlet piece. Consider a specihc gg-pair in |c) that is now in 
a common representation, for example i and j in eq. (A.12) being in the representation a. If the qq-paii 
is also paired into a common representation in the basis vector, the quark trace gives a factor of Tji, Nc 
or vanishes. The resulting color structure is topologically equivalent to a color structure of the same order 
as eq. (A.9) with i and j exchanged for a gluon or singlet, apart from possibly being disconnected (due to 
the hrst term in eq. (A.ll)) and containing vertices (from the second term of eq. (A.ll)). Such color 
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structures are then equivalent to color structures with one less (;g-pair, and hence the (/g-pair ij does not 
require any other manipulation than what has already been shown. However, if i and j are not paired in 
the basis vector, i.e., they do not belong to the same qq pair of the basis vector, the color structure is not 
guaranteed to contain loops of the type in eq. (A.l) or eq. (A.7). Hence additional steps are required. 

The terms in eq. (A. 13) when contracted with the basis vectors will be of the form 



(A. 14) 


where the gray blob contains everything else in the vacuum bubble (the rest of the chain of arbitrary 
representations and the blob in eq. (A.13)) and a, b and c label the representations of the thin double lines, 
singlets or octets. The goal is to manipulate this color structure into a color structure where the quark-line 
from / to /c is only connected to the basis vector by one representation, a singlet or an octet. If this is 
achieved, the quark trace that I and k are part of contains fewer generators, and the procedure can be 
repeated until the trace is over three or fewer generators. There are two possibilities for c in eq. (A. 14): 

1. If c is a singlet: 

(a) and at least one of a and 6 is a singlet, then I and k can be seen as a quark-antiquark pair that is 
paired in the basis vector. The quark-antiquark pair Ik is in a singlet if both a and h are singlets, 
and otherwise in an octet. 

(b) and both a and h are octets, then eq. (A.11) can be used to combine a and 6 in a singlet or an 
octet. When a and h are in a singlet they are part of a loop of the form of eq. (A.7), and when 
they are in an octet they form a loop as in eq. (A.l). Removing these loops will leave a color 
structure where I and k are in a singlet in the first case, and in an octet in the second case. 

For both (a) and (b), the vacuum bubble is smaller, since the basis vector contains one quark-antiquark 
pair less; initially it had the pairs li and jk, but they have been replaced by one quark-antiquark pair 
Ik. 

With Ik being in a singlet or an octet and connected to the basis vector, the procedure can be repeated. 
The color structure can be written on the form of eq. (A. 14), with Ik taking the place of either li or 
jk, until the new c is an octet or the new I and k are directly connected, in which case the quark 
trace is over two, or fewer, generators. 
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2. If c is an octet the color structure is of the form 



(A.15) 


We wish to manipulate the color structure such that, instead of a, b and c, there is only one repre¬ 
sentation connecting to Ik and it connects directly to the basis vector. 

Using a completeness relation, eq. (3.2), on quark I and antiquark j in eq. (A.15) gives 



where the first term vanishes if a is a singlet (since the generators are traceless). For the first term, 
there is only one representation, b, connecting to the quark-line Ik, so the procedure can be repeated 
to remove the rest of the octets/singlets from the trace. For the second term of eq. (A. 16) another 
completeness relation can be applied to quark I and antiquark k, 



(A.17) 

The first term of eq. (A.17) vanishes unless b is an octet. As the quark trace with I and k has been 
disconnected they can be seen as a singlet connected to the chain of arbitrary representations and the 
procedure can be repeated. For the second term, the different possible combinations of a and b can 
be divided into two cases. 
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(i) The representations a and b are octets. The traces over three generators can be replaced by if 
and d vertices and a loop of type eq. (A.5) is formed 



(A. 18) 


where the two small gray blobs are different combinations of if and d, and g labels the gluon 
connecting to I and k. Contracting the loop of the type in eq. (A.5) leaves a color structure 
where g connects the basis vector to I and k, and c connects the big gray blob to the basis vector. 
The color structure is now topologically equivalent to a color structure with one less g'^-pairs, 
and hence the procedure can be repeated. 

(ii) For the case when a is a singlet (analogous steps can be performed if 6 is a singlet) and b is 
either a singlet or an octet, a completeness relation, eq. (3.2), can be applied 



E 



If 6 is a singlet, then two gluons have been connected to the chain of arbitrary representations, 
but as a and b are both singlets, there are now as many gluons connected to the chain of arbitrary 
representations as if a and b had been octets instead. This is important for the latter part of 
this appendix, where constraints are put on the representations in the Wigner 6j coefficients. 
Note that if both a and b are singlets, then the completeness relation can be applied to any 
representation in the chain of arbitrary representations, not just /3. If b was an octet, the trace 
over three generators can be written as if and d vertices and the resulting vertex correction can 
be removed as in eq. (A.3). 

In both cases, (i) and (ii), the quark trace is shorter, and the procedure can be repeated. 

All possible cases have been covered above, the quark traces are systematically removed by disconnect¬ 
ing representations from them, making them shorter until they are over three, or fewer, generators. After 
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the quark traces have been removed, the color structure can be contracted using the types of loops required 
for a gluon-only amplitude of the same order as the original amplitude with quarks (with the addition 
of loops of type eq. (A. 7)). The procedure in this section of the appendix is valid for arbitrary order in 
perturbation theory, in particular for LO and NLO. By allowing contraction of loops of the type in eq. (A. 5) 
any LO and NLO QCD color structure can be completely contracted. 

A.3 Constraints on representations in Wigner coefficients 

The above reasoning proves that only Wigner coefficients of the form of eq. (4.2) and eq. (4.3) are required, 
since those are the only forms occurring in eq. (A. 2 ), eq. (A. 3), eq. (A. 6 ) and eq. (A. 8 ). We will now move 
on to put constraints on the representations appearing for LO gluon-only color structures, i.e., in eq. (A. 2), 
eq. (A.3) and eq. (A. 8 ). After this, constraints on the representations required for up to NLO with both 
gluons and quarks are handled, i.e., eq. (A. 6 ). 

An important property of the representations in tensor products between multiple adjoint representa¬ 
tions is their so-called first occurrence [20]. The first occurrence, nj, of a representation M is the lowest 
integer i such that M G A®*. In [20] it was shown that uj of the representations \n M ® A can differ from 
Uf of M by at most ±1. This is a key part in the following derivation of constraints on the representations 
in the Wigner coefficients. 

We want to show that the above contraction procedure never requires tp in eq. (A.2) to have a first 
occurrence larger than [n/2j for tree-level color structures. To prove this, we note that it is enough to 
prove it for the basis vectors with highest possible first occurrence for all involved representations. The 
representations in the basis vector are 

01,02, • • • ,a[n/2j-2,«LV2j-l>“Ln/2j, • ■ ■ , On-S 
and the highest first occurrences they can have are given by 

2.3.. [|J - 1, [|J , [|J - neveii 

2.3.. ..,[5J-l,[eJ,[iJ,[eJ-l,...,2 nodd. 

The rightmost color structure of eq. (A.2) will be zero if nf{ip) > u/(o) -|- 1 or nf{^jJ) > nf{'y) -|- 1, since 
the representations can differ by at most ±1 in first occurrence. From these limits we see that the first 
occurrence of ip can only exceed [n/2j if nf{a) > [n/2j and nf{'y) > [u,/2j. Noting that in the chain of 
representations o and 7 are not next to each other (as they are separated by /3), one can see from eq. (A.21) 
that o and 7 cannot both have first occurrence [n/2j since they are not adjacent. Thus the first occurrence 
of Ip cannot exceed [n/ 2 j. 

Stronger constraints can be put on the coefficients by examining the first occurrences of a, /?, 7 and ip 
of eq. (A. 2). For even n, we see from eq. (A.21) and eq. (A.2) that the representation ip can only have first 
occurrence of [n/2j if nf{a) = nf{'y) = [n/2j — 1. The consequence of this is that for even n there can be at 
most two representations in coefficients of the form of eq. (4.3) with first occurrence [n/2j, and they cannot 


(A.20) 


(A.21) 
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be adjacent in the coefficient (in eq. (A.2) they would be (5 and ip, thus separated by a and 7 ). If n is odd, 
then two of a, /? and 7 can have first occurrence [n/ 2 j, and ip may also have first occurrence [n,/ 2 j, the 
Wigner coefficient can then contain at most three out of four representations with first occurrence [n./2j. 

For the Wigner coefficients of the form of eq. (4.2) in eq. (A.3) the first occurrences are the same as 
for three adjacent representations in the chain of arbitrary representations, eq. (A.20), i.e. for even n there 
can be at most one representation with [n/2j and for odd n there can be at most two. For the Wigner 
6 j coefficient occurring in eq. (A. 8 ), the same argument applies, no new representation has been added so 
it cannot contribute to an increase in first occurrences of representations. These constraints were used in 
calculating the number of required Wigner coefficients for LO gluon-only processes, shown in table 2. 

It is now of interest to see how the constraints on the representations in the required Wigner coefficients 
change by including NLO diagrams and quarks. As the only change compared to the previous discussion 
is eq. (A. 6 ), it is only the constraints on the Wigner 6 j coefficients with three general representations, 
eq. (4.2), that are altered. The case of interest is when at least two of the three representations have a first 
occurrence of [n/ 2 j (the case with one representation with Uf = [n/ 2 j is already required for LO gluon 
color structures), this case occurs when ip has a first occurrence of [n/2j. For even n, the highest first 
occurrences appear when nf{a) = nf{'y) = [n/ 2 j — 1 and nf{f3) = [n/ 2 j, for which we see from the right 
hand side of eq. (A. 6 ) that coefficients with two representations with first occurrence [n/ 2 j may appear {ip 
and /3). Similarly, for odd n, the interesting case is when nf{a) = nf{l3) = [n/2j or nf{l3) = nf{^) = [n/2j 
for which all three representations can have first occurrence of [n./2j. As is shown in table 2, these altered 
constraints lead to a very small difference in the required number of coefficients between tree-level gluon 
color structures and NLO color structures with quarks and gluons. 


B Construction history independence of Qj coefficients 


To prove that the 6j coefficients which only contain unique vertices are independent of the construction 
history, we note that the various copies of a given vector space, such as the and the C^^^’^^-versions of 

the vector space corresponding to the 35-plets for three gluons, are isomorphic. With unique vertices 
we mean vertices between representations q;,/ 3 and 7 , s.t. the tensor product a 0 /3 only contains one 
instance of the representation 7 . 

Since the vector spaces 1 /^’^’^® and l/io,35 isomorphic, we can find a transformation U from to 
A®^, which maps the vector space 1 /^^’^® to ^ unitary way, but which maps every other irreducible 

subspace in A®^ to 0. Seen as a matrix (in a suitable basis), U thus has an off-diagonal 35-dimensional 
matrix block whereas all other elements equal zero. The projection operators pio,35 p27,35 gj.g then 
related as _ uP^^-SS-y-i^ birdtrack notation. 


-p 10.35 






TT-l- 



u 




u 


(B.l) 


For the proof of the existence of construction history independent 6j coefficients, the existence of U is 
enough, but for explicit construction, we remark that U can be obtained by noting that any invariant 
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transformation T between the irreducible vector spaces yio,35 y27,35 |^g proportional to the unit 

matrix by Schur’s lemma. Since any transformation expressed in birdtracks is invariant, we can use any 
non-vanishing T. Normalizing T and inserting it in the only non-vanishing block in U gives us U. 

Using U we can rewrite the 6j coefficient 




from which we see that if the vertex 7;">/^'(27,35) jg (defined s.t. it equals the encircled region in 



^ ^q,/3,(27,35)^ 


(B.2) 


(B.3) 


and similarly for the vertex to the left in the 6j coefficient, the equality is satisfied. For unique and 
normalized vertices, the r;"’^’(^'^’^^)-vertex has by necessity already been defined this way (up to a sign), 
and we conclude that the 6j coefficients must agree. By the same argument all 6j coefficients involving 
unique vertices which only differ by the construction history must agree (modulo signs). For vertices 
appearing in several forms, or to fix the sign ambiguity, we note that if we define the various vertices s.t. 
eq. (B.3) holds, the 6j coefficients agree. 


C Proof of symmetries of 6j coefficients 


The rotation symmetry of the first type of coefficient in eq. (4.9) is obvious. For the second type of 
coefficient, the Wigner coefficient can be put into a form where the symmetry is easier to see. Moving the 
central vertex upwards and to the right, rotating and finally moving the upper right vertex to the middle 
gives the symmetry. 





a 


(C.l) 











































To prove the conjugation symmetry, eq. (4.10), we recall that the color sum in the Wigner coefficients 
only consists of factors with a different number of closed quark-lines, i.e., they can be expanded in the 
fundamental representation, giving a polynomial in Nc (up to vertex normalizations). This means that the 
Wigner coefficients are real numbers, and as such do not change under conjugation. 

The final symmetry, eq. (4.12), is shown by simply swapping places of the lower left and lower right 
vertices. 
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